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Abstract 



Determination of the orbital angular momentum of the proton is a difficult but 
, important part of understanding fundamental structure. Insight can be gained from 

Q I suitable models of the gluon asymmetry applied to the = 1/2 sum rule. We have 

constrained the models of the asymmetry to gain possible scenarios for the angular 
momentum of the protons constituents. Results and phenomenology for determining 
Lz are presented. 



1 Status of Proton Spin Structure 

For the past twenty years, much work has been done to understand the spin structure 
of the nucleons. There has been progress in determining the contribution of the lightest 
quarks to the spin, but there is still uncertain knowledge about the gluon contribution. 
Transversity studies have contributed additional insight about quark dynamics, but little 
is known about the the orbital angular momentum of the constituents. [1] This paper will 
summarize a project that provides a method of gaining insight into the nature of the 
orbital angular momentum of the nucleon constituents. 
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Recent experiments [2, 3] have significantly lowered the measurement errors of the 
quark longitudinal spin contribution (AE) to the proton. The COMPASS collaboration 
analysis quotes a result 

AE = 0.30±0.01(stat) ±0.02(evol), all data (1) 

while the HERMES collaboration analysis quotes a result 

AE = 0.330 ± 0.025(exp) ± 0.011 (th) ± 0.028(evol), all data. (2) 

These groups and others [4] have been working on providing a significant measure of the 
proton's spin weighted gluon density, 

AG{x,t)=G++{x,t)-G+4x,t), (3) 

where x is the Bjorken scaling variable and t = log{as{Ql) I loglo'slQ^)) is the evolution 
variable. The combination of these measurements is summarized in terms of the Jz = \ 
sum rule: 

J, = ^ = ^AE + AG' + L,. (4) 

Here AE = /q^ dxAqix^ t) and AG = Jq dxAG{x, t) are the projections of the spin carried 
by all quarks and the gluons on the z-axis, respectively. Also is the net z-component 
of the orbital angular momentum of the constituents. We do not attempt to separate the 
fiavor components of within the sum rule. 



2 Modeling the Gluon Asymmetry 

Experimental groups at the COMPASS, HERMES and RHIC collaborations are measur- 
ing both the gluon polarization and the asymmetry, A = AG/G to determine the gluon 
polarization [2, 3, 4]. Since there is no suitable theoretical model for AG, we have devised 
a way to model the asymmetry, A{x, t) to gain insight into the structure of AG. This, 
coupled with the Jz = \ sum rule can then shed light on the nature of the orbital an- 
gular momentum of the constituents, L^. To model A{x,t), we write the polarized gluon 
asymmetry using the decomposition 

A{x,t) = AG/G = Ao{x) + e{x,t), (5) 

where 

M.) ^ [(-^)/( 3^)] (6) 

is a scale invariant calculable reference form [5]. Here e{x,t) represents the difference 
between the calculated and gauge-invariant asymmetry. Since AG is unknown, a useful 
form is to write equation (5) as 

AG = Ao{x) G{x, t) + AGe{x). (7) 

Although the quantity AGe{x) is not a physical parameter, it allows the theoretical devel- 
opment of the calculable quantity, Aq. Once an asymmetry is generated from equations 
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(6) and (7), the gauge-invariant quantity A{x,t) can be compared to data. Thus, each 
Ansatz for AG^{x) gives a corresponding form for AG and a parametrization for L^. 
These can be compared to existing data to provide a range of suitable models for these 
contributions. 

With the definition for the asymmetry in equation (6), the DGLAP equations can 
then be used to evaluate the evolution terms on the right side. 



^0 



APcq O Ag + APgg ^ AG 
PGq^q + Pgg ® G 



(8) 



The polarized gluon distribution in the numerator of equation (8) is replaced by AG = Aq- 
G+ AGe. For certain unpolarized distributions, there are points at which the denominator 
vanishes. To avoid this, we write equation (8) as: 

dAG 



dt 



= (2//5o) 
= Ao 



AP, 



LO 
91 



(Ao-G + AG,) 



(9) 



dC 
'dt 

= (2//3o)^ 



oLO 
91 



G 



The NLO form is essentially the same as equation (9) with the splitting functions P^^ 
replaced with their NLO counterparts. The quark and gluon unpolarized distributions 
are CTEQ5 and the polarized quark distributions are a modified GGR set. [6] 

There are constraints on Aq{x) that must be imposed to satisfy the physical behavior 
of the gluon asymmetry, A[x). These are: 

• positivity: |>lo(2^)| < 1 for all x, and 

• endpoint values: Aq{{)) — {) and Aq{1) — 1 

Note that the constraint oi Aq ^ 1 is built in to satisfy the assumption that the large x 
parton distributions are dominated by the valence up quarks in the proton. The convo- 
lutions are dominated by the quark terms, which force the asymmetry to unity as a; — > 1. 
To investigate the possible asymmetry models, we use a parameterization for Aq in the 
form 



Ao = Ax"" -{B-l\ 



X' 



/3 



[B - A)x 



/3+1 



(10) 

which automatically satisfies the constraints that Ao{0) = and Ao{l) = 1. Once a 
parametrization for AGe{x) is chosen, equation (9) is used to determine the parameters 
in equation (10). 



3 Results and Conclusions 

The models for AG^Xx) that led to asymmetries that satisfied these constraints were all in 
the range | AGedx\ < 0.25, with positive and negative values included. Larger values 
of AGe violate one or both of the constraints. A representative sample of models that 
satisfy the constraints are hsted in Table 1. 

Note that the integrals for AG are all positive, ranging from about 0.01 to 0.42. 
The models that gave negative values for these integrals did not agree with the existing 
asymmetry data, reported at this workshop to be: 
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Table 1: Gluon Asymmetry Parameters 







A 


Jg AAtj^UX 








o 1 "i o 2 2 1 S 2 


0.18 


2(1 — x) 


0.25 


/(16 A 21 t 31 

4x — 4x + X 


0.42 


-2{i-xy 


-0.25 


l.lbx^-^ - 1.5x^-1 + 0.75x^-1 


0.01 


-90x^(1 - xY 


-0.25 


3.5x^-3 - 4.5a;2-2 + 2a;3-2 


0.05 


9x{l - xY 


0.125 


3.75x1-^ - 3x1-6 + 0.25x2-6 


0.29 


-9,x(l -x)'^ 


-0.125 


3.25x^-4 - 3.75x2-2 ^ ^ 5^3.2 


0.11 


4.5x(l - xY 


0.0625 


2a;°-9 - 1.5a;i-2 + 0.5a;2-2 


0.37 


-A.hx{l-xY 


-0.0625 


2.25x1-1 - 2.25x1-9 + x2-9 


0.23 



• AG/G = 0.016 ± 0.058 ± 0.055 at x = 0.09 from COMPASS, > 1 GeV^ 

• AG/G = 0.060 ± 0.31 ± 0.06 at x = 0.13 from COMPASS, < 1 GeV^ 

• AG/G = 0.078 ± 0.034 ± 0.011 at x = 0.204 from HERMES, factorization method 

• AG/G = 0.071 ± 0.034 ± 0.010 at x = 0.222 from HERMES, approximate method. 

The models in Table 1 that are within one a of the preliminary data stated above are in 
the third, fourth and sixth ro-ws, respectively. Plots of the full asymmetry are shown in 
Figure 1. None of the models in Table 1 are ruled out by the data since they fall within 
two a of the data for our values of > 1 GcV^. All of these models except for the 
fourth row in Table 1 (impulses in figure 1) generate total asymmetries A{x,t = 0) that 
are close to A{x) = x. Ironically, early assumptions of the polarized gluon assumed this 
functional form as a naive estimate to the asymmetry. Next-to-leading order corrections 
to these asymmetries tend to bring them less positive, but with the same general shape. 
A full set of viable asymmetries will be presented in an upcoming paper. [7] 

Using the data on AE in Section 1, the relation between < AG > and < > can be 
written as: 

< AG 0.35 - <L^> ±0.02. (11) 

The three models of the asymmetry that agree most closely with existing data give values 
of AG in the approximate range of — > 0.11. Thus, the existing data with equation 
11 imply the approximate relation 0.24 < < 0.35 ± 0.02. Thus, the contribution of 
the orbital motion of the constituents to the proton spin may be comparable to the total 
quark contribution. A recent lattice calculation of the contribution of the quark orbital 
motion to the proton spin [Lf) is consistent with zero. [8] Thus, the gluonic orbital 
motion appears to provide the majority contribution to in the Jz — \ sum rule. It is 
clear that future measurements of AG and AG/G must be made in a wider kinematic 
range of x and with improved precision to better specify the appropriate model of the 
asymmetry and to extract the x and dependence of the orbital angular momentum of 
the constituents. 
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Figure 1: The gluon asymmetries most closely in agreement with data. Solid line, impulses 
and linespoints represent the models in rows 3, 4 and 6 of Table 1 respectively. 
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